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Abstract. We discuss the constraints on the AGN evolu- 
tion from the cosmic X-ray background and source counts. 
A synthesis model to fit the X-ray background is pre- 
sented. The data used are a compilation from Gruber 
(1992), mainly based on HEAO-1 data. The spectrum of 
type 2 AGN has been modeled including Compton down- 
scattering within the absorbing material. We found an im- 
provement when a dependence on redshift of the relative 
number of obscured sources is introduced. In particular, 
we found a decrease of the fraction of type 2 AGN at red- 
shifts larger than ~2. 

However, only the soft X-ray counts are well repro- 
duced, while the hard (5-10 keV) BeppoSAX/HELLAS 
source counts are underestimated by a factor of ^^1.5. 
A satisfactory global solution, including the good re- 
production of the hard counts, is instead obtained if 
the XRB normalization, recently obtained from Bep- 
poSAX/MECS measurements (about 30% higher than 
HEAO-1), is adopted. 

Key words: X-rays: general - galaxies: nuclei - luminos- 
ity function - cosmology: diffuse radiation 



1. Introduction 

The cosmic X-ray background (XRB) above ^ 1 kcV 
is known to be produced by integrated emission of dis- 
crete sources. XRB synthesis models are usually based on 
the so-called unification scheme for Active Galactic Nuclei 
(AGN) , which ascribes the different observational appear- 
ances to the orientation of a molecular torus surrounding 
the nucleus. The intersection of the line of sight with the 
torus determines a type 2 AGN (AGN2), while the direct 
observation of the nucleus identifies a type 1 AGN (AGNl) 
(Antonucci & Miller 1985). An appropriate combination 
of the two type of sources can reproduce the main char- 
acteristics of high-energy XRB, i.e. intensity and spectral 
shape (Setti & Woltjer 1989; Madau et al. 1994; Matt & 
Fabian 1994; Celotti et al. 1995; Comastri et al. 1995). 



In this paper, the constraints on the evolution of AGN 
from the XRB and the source counts are examined, and 
a possible variation of the standard model (Comastri et 
al. 1995), in which the AGN2/AGN1 ratio changes with 
redshift, is explored. We made use of source counts data 
in the 0.5-2 keV band by ROSAT (Georgantopoulos et al. 
1996), in the 0.3-3.5 keV band by Emstein (EMSS; Gioia 
et al. 1990), and the source counts provided by the Bep- 
poSAX HELLAS survey in the 5-10 keV band (Fiore et 
al. 1999; Comastri et al. 1999). The XRB fitting proce- 
dure was based on the Gruber (1992) compilation of data, 
mostly from HEAO-1 in the 3-50 keV range (Marshall et 
al. 1980). 

Throughout this paper, a Hubble constant Hq — 50 
Km s~^ Mpc~^ and deceleration parameter qq — 0.5 have 
been assumed. 

2. Models for X-ray spectra and evolution of AGN 

2.1. AGNl and AGN 2 spectra 

The local AGN spectrum is assumed to be the sum of 
AGNl and AGN2 spectra {Fi{E), F2{E) respectively) 
weighted by the number ratio of AGN2 to AGNl, R: 

FiociE) [Fi{E) + R X F2{E)] . 

Following Comastri et al. (1995), a double power-law 
with a Compton reflection component Fr{E) has been 
adopted for the AGNl spectrum: 



Fi{E) oc 



E < 1.5 keV 
Fr{E) E>1.^ keV 



with a2=0.9 (Matsuoka et al. 1990; Pounds et al. 1990; 
Nandra & Pounds 1994), and ai=1.3. 

The steeper, low-energy {E < 1.5 keV) spectrum rep- 
resents the so-called soft excess. The shape and contribu- 
tion of this component is not well known, and in many 
sources evidence for its very existence is lacking alto- 
gether. As a baseline, we have adopted the same prescrip- 
tion as Comastri et al. (1995), but as it now appears to be 
rather extreme, we have also explored the opposite case, 
i.e. no soft excess at all (see Sect. 3.1). 
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Fig. 1. The AGNl spectrum (solid line) and the AGN2 average 
spectrum (dashed line) as produced by the model. 



The adopted value of the cut-off energy Ec=400 keV 
is also similar to those used in previous models (Comastri 
et al. 1995; Celotti ct al. 1995) even if recent BeppoSAX 
results (Matt et al. 1999a and references therein) seem 
to suggest somewhat smaller values, but with a rather 
large spread. The number of sources with reliably mea- 
sured values is so low, however, that we preferred to still 
use the values adopted in previous models, to make easier 
the comparison. In any case, we tested the effect of adopt- 
ing lower values for the cut-off, and found that the peak 
in the XRB spectrum at ~ 30-40 keV is less well fitted. 

The term Fr{E) represents the Compton reflection 
component by the accretion disk and by the torus in- 
ner surface and has been evaluated following Magdziarz 
& Zdziarski (1995), assuming an inclination angle of 60°. 

According to unified schemes, AGN2 spectra arc ob- 
tained as AGNl spectra seen through absorbing mat- 
ter. The distribution of equivalent hydrogen column den- 
sity {Nh) is chosen to be logarithmic, i.e. '^^[^"^^^^ oc 
log Nh, that is a reasonable analytical approximation to 
the recent data on Seyfert galaxies (Maiolino et al. 1998; 
Risaliti et al. 1999). 

As described in a previous paper (Matt et al. 1999b), 
we developed a transmitted spectrum model by means of 
Monte Carlo simulations, assuming a spherical geometry 
with the X-ray source in the centre and considering photo- 
electric absorption, Compton scattering and fluorescence 
(for iron atoms only), fixing element abundances as tabu- 
lated in Morrison & McCammon (1983). 



This transmitted component, which is relevant for 
Nh > lO^'^cm^^, has been so far included in XRB synthe- 
sis models only by a handful of authors (Madau et al. 1994; 
Celotti et al. 1995; Matt et al. 1999b; Wilman & Fabian 
1999). The final spectrum Fn^{E) has been then aver- 
aged over the iV//-distribution to obtain the total AGN2 
spectrum: 

^ ^ I!:^n::PnAE) log Nh dHog Nh) 
C log Nh dilog Nh) 

where the TV^f-distribution has been considered in the 
range lO^^ < Nh < 10^^ cm-^. 

The AGNl and AGN2 spectra are shown in Fig. 1. 
The AGNl spectrum is flattened by the reflection compo- 
nent, whose contribution reaches ~ 1/3 of the total at its 
maximum {E ~ 30 keV). The addition of a transmitted 
component significantly affects the total AGN2 spectrum, 
again with a flattening around 30-40 keV. The Compton 
scattered photons in AGN2 increases the total spectrum 
by ^ 20% at ~ 30-40 keV with respect to a model involv- 
ing only absorption. 

The last step in evaluating the overall local spectrum 
concerns the choice of the number ratio R. Maiolino & 
Rieke (1995) find i? = 4.0 ± 0.9, if type 1.8, 1.9 and type 
1.2,1.5 Seyfert galaxies are respectively classified as AGN2 
and AGNl: the estimate agrees with previous (Huchra & 
Burg 1992; Goodrich et al. 1994) and more recent (Ho et 
al. 1997) results, and so R = 4.0 has been adopted. 

For the sake of simplicity, iron emission line has not 
been included, even if it is a common feature in AGN. 
However, the contribution of the line to the 1.5-7 keV 
XRB is expected to be less then 7% (Gilh et al. 1999a) 
and it is smeared out by the integration over the redshift 
range, so that XRB retains its characteristic smoothness 
(Schwartz 1992), unless the emission is dominated by a 
small range of redshifts (Matt & Fabian 1994). 

2.2. Cosmological evolution 

Hard X-rays (3 > keV) are well suited for the selection 
of type 1 and, especially, type 2 AGNs as they are less 
affected by absorption. Until recent past little was known 
about the evolution of the AGNs in this band. The data 
obtained from the ASCA satellite have allowed a first de- 
termination of the 2-10 keV AGN XLF (Boyle et al 1997). 
However the statistics was still poor. On the contrary the 
AGNl XLF in the soft X-ray band is retained to be well- 
known at low and intermediate redshift (Boyle et al. 1994; 
Page et al. 1996; Jones et al. 1997), while at higher red- 
shift {z > 3) insufficient sampling and lack of statistics 
prevent the XLF from being firmly evaluated. Anyhow, 
X-ray AGNl are detected up to z = 4.6 (Miyaji et al. 
1998) and no evidence for a space density turn-over is 
found up to ~ 3, likewise in the optical (e.g. Kennefick 
et al. 1996) and radio surveys (Shaver et al. 1999). 
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For these reasons we chose to tie the AGN evolution 

to the soft X-ray XLF of AGNl. Wc used the pure himi- 
nosity evolution (PLE) scenario which fits the combined 
ROSAT and EMSS data on AGNl space density (Boyle et 
al. 1994). This corresponds to a local luminosity function 
that can be represented by a double power-law where the 
break-luminosity , i.e. the luminosity value correspond- 
ing to the slope change, evolves as L,(z) oc {1 + zy. In 
the following, we adopt the PLE H-model of Boyle et al. 
(1994) (hereafter B94) in the 0.3-3.5 keV band: 

with 71 = 1.36, 72 = 3.37, = 0.39, = 1.45 x lO"'' 
Mpc-3 (10"^^ erg s'^)-^ and L44, in unity of 10^^ erg 
s The break- luminosity evolution follows: 

^ = 0) X (1 -h^)*" ^ Z < Zmax 

\L^[z = 0) X (1 -|- Zmax) Z > Zmax 

where k = 2.90 and Zmax = 1-73. We limited the anal- 
ysis of AGNl XLF to a PLE model because the most 
recent attempts with pure density evolution models, in 
which the space density $(i) directly evolves in redshift, 
overproduce the soft XRB (Hasinger 1998). We have also 
to introduce an AGN2 XLF, which is a matter of strong 
debate. In the framework of the unification scheme, we 
assumed the density of AGN2 to be i?=4 times that of 
the corresponding unobscured AGNl. We assumed the 
A''ij-distribution to be independent of the AGNl source 
luminosity. A different approach may consist in setting a 
completely unrelated XLF, but it would involve too many 
parameters and there are not enough data to yield a reli- 
able estimate. 

Boyle et al. (1997) directly measured the AGNl and 
AGN2 XLF on a sample of 26 2-10 keV ASCA sources at 
a flux limit of 5 x 10~^^ erg cm~^ s~^, combined with the 
HEAO-1 AGN. The analysis gave a result consistent with 
the 0.3-3.5 keV AGNl XLF, albeit the evolution seems to 
be slower {k = 2.04). 

3. The XRB synthesis model 

The intensity of the XRB has been calculated as follows: 

ci^N F[L,E{l + z)] 



m 



where 



dVdL 



Andf 



'-dLdV 



dVdL 



is the AGNl XLF, di the luminosity dis- 
tance, F[L,E{1 -\- z)] the AGN spectrum. Introducing 
the comoving volume dV{z) and the spectrum normal- 
ization K corresponding to the considered luminosity, i.e. 
K J^j^Fi{E)dE = L (A£;=0.3-3.5 keV), the above rela- 
tion can be written as: 



I{E) 



Kc 



d^N Fio,[E{l + z)] 

'«^^c;^(i + 0)'(i + 2go^)5 



where the integration is performed in the range 10 < 
L < lO^'^erg s-\ according to the B94 0.3-3.5 keV XLF. 

The XRB synthesis has been used to explore the evo- 
lution of the AGNs. The model was fitted to the data 
discussed by Gruber (1992), largely based on HEAO-1 A2 
3-50 keV data^. 

Best fit parameters were determined by a x^- 
minimization procedure based on the CERN MINUIT 
software package. 



3.1. The baseline model 

The first analysis concerns the AGNl high-redshift XLF. 
We tried to model the decrease of AGNl at high redshift 
by introducing a new parameter Zd > Zmax corresponding 
to the onset of the space density decrease. For redshifts in 
the range Zmax < z < Zd, the density stays constant. We 
used three different shapes for z > z^: i) an exponential 
{N{z) oc e-"(^-^<i)); 11) a polinomial {N{z) oc (^)-"); 
III) a sharp cut-off {N{z) = 0). We allowed a to vary, and 
fixed Zd=4:.5 (which is a simple but still good representa- 
tion of the space density of AGNl according to the present, 
sparse observations). For any of the laws described above, 
the XRB is overproduced. We then allowed also Zd to vary, 
always obtaining values of Zd < 3. This is a value much 
lower than observed as discussed in Sect. 2. 2, and therefore 
all these models can be considered unacceptable. 

We thus decided to adopt a sharp cut-off (model ill 
above) with Zd = 4.5 and to fit the XRB changing a dif- 
ferent parameter. We first tried to vary the number ratio 
R, obtaining a value of R = 3.5 ± 0.1; the error corre- 
sponds to 68% confidence level, following Lampton et al. 
(1976). The value is still in agreement with the observa- 
tions. The fit was not very satisfactory from a statistical 



point of view, = 



41 



Part of the high value of the 



is certainly due to the fact that we have not included other 
classes of sources, notably starburst galaxies and clusters 
of galaxies, which are likely to contribute at low energies. 
However, the fit is bad also at energies where we expect 
the AGN to dominate. The fit shows also an excess at 
higher energies {E > 70 keV), probably due to the value 
assumed for the intrinsic cut-off energy; however, as un- 
derlined in section 2.1, the choice was motivated by the 
aim of reproducing well the ~ 30-40 keV peak. Even if we 
know that a really good fit is impossible with any smooth 
model, because of the significantly large fluctuations of 
the data (mostly due to the fact that measurements from 
different instruments have been used simultaneously), we 
decided to try to improve the quality of the flt by further 
changing our baseline model. 



-dLdz 



^ Recently, a new XRB data set by the UCSD/MIT Hard 
X-Ray and Gamma- Ray instrument (HEAO-1 A4) has been 
published (Gruber et al. 1999). These data join smoothly with 
the XRB data reported by the GSFC HEAO-1 A2 instrument, 
used in this paper. 
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3.2. The R(z) model 
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Fig. 2. The number ratio as a function of redshift R{z) in 
unity of the local value Ro, fixing Ro = 4.0 (solid line). The 
dashed lines enclose the 68% confidence region. 



• : Gruber (1992) data 


X : ASCA data 
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Fig. 3. The fit to the XRB data (solid Une) with the R{z) 
model; ASCA data are also reported, even though they have 
not been considered in the fit. 



Because the evolution of type 2 AGNs is probably the 
least known ingredient of the model, we tried to intro- 
duce a z-dependence of the AGN2/AGN1 ratio, R. Dif- 
ferent analytical models have been tried for R{z): mono- 
tonic (both exponential and power law) shapes and King 

profile (i?(z) oc ^1 -|- ), but none of them gave a 

less than 

We found, instead, a better result when a two param- 
eter function of the form: 



R{z) = i?o X (1 z) 



is used, with a significant improvement of the fit (x^ = 
2j^) (see Figs. 2 and 3). According to the F-test, this 
corresponds to a 95% confidence level with respect to the 
i?o=constant fit described above. The best-fit parameters 
are: ai = 1.28 ± 0.04, a2 = -0.89 ± 0.03. The decrease of 
the x^ is mainly due to a better reproduction of the data 
in the 5-10 keV energy range and around the XRB energy 
density peak (~30 keV). 

We then checked how much dependent the above re- 
sult is on the adopted shape of the local AGN spectrum, 
and in particular on our choice for the soft excess which, 
as discussed in Sect. 2.1, may be rather extreme. To do 
that, we adopted the opposite assumption, i.e. absence of 
a soft excess. The result gives an even quicker decrease of 
AGN2 with redshift. This is due to the fact that in ab- 
sence of the soft excess the hard power law has an higher 
normalization in order to produce the same 0.5-3.5 keV 
luminosity, which implies higher luminosity at harder en- 
ergies. Thus, in order not to overproduce the XRB, less 
AGN2 are required. 

The effect of using a pure photoelectric absorption 
model for AGN2, as done in some previous papers, has 
also been tested. In this case the spectrum intensity is 
higher than the photoelectric absorption -|- Compton scat- 
tering model (Matt et al. 1999b), and again a steeper R{z) 
emerges. 

It is important to note that the introduction of R{z) 
implies the presence of a density evolution component for 
AGN2. Density evolution has been taken into account by 
the ROSAT Ah Sky Survey (RASS) 0.5-2 keV data analy- 
sis of Miyaji et al. (1998), who introduce a luminosity de- 
pendent density evolution (LDDE) basically correspond- 
ing to a sudden drop in the evolution rate at lower lumi- 
nosity. However, the RASS sample contains both AGNl 
and AGN2, so the authors cannot outline any differences 
between the two kind of sources and estimate their evolu- 
tionary properties separately. 

So far we derived the AGN evolution only by fitting 
the XRB. Then we checked the consistency of our re- 
sults with the soft (ROSAT and Einstein, 0.3-3.5 keV) and 
hard (BeppoSAX, 5-10 keV) counts. In Fig. 4 the source 
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Fig. 4. (Upper panel) The integral soft counts (0.3-3.5 keV) 
compared to ROSAT -|- EMSS data. (Lower panel) The integral 
hard counts (5-10 keV) compared to HELLAS data. The source 
counts have been evaluated including the R{z) term. 



counts predicted by the model in the two bands are com- 
pared to the ROSAT (Georgantopoulos et al. 1996) -|- Ein- 
stein/EMSS (Gioia et al. 1990) and BeppoS AX/HELL AS 
(Flore et al 1999) data. 

A good agreement in the soft band is found, while the 
hard counts are underestimated by a factor of ~ 1.5 at ~ 
10~^^ erg cm~^ s~^. A possible solution of this discrep- 
ancy is discussed in the next section. 



Fig. 5. Same as in Fig.2, but introducing a 30% higher nor- 
malization to XRB data. 

3.3. The normalization 

Recently, new XRB data up to 10 keV have been obtained 
with imaging instruments on-board ASCA (Gcndrcau et 
al. 1995, Miyaji et al. 1998) and BeppoSAX (Molendi ct 
al. 1997, Parmar et al. 1999). While there is still some dis- 
agreement between the XRB normalizations obtained by 
BeppoSAX and ASCA (Parmar et al. 1999), both of them 
arc higher than that obtained by HEAO-1. For instance, 
BeppoSAX/MECS result is 30% higher than the HEAO-1 
(Vecchi et al. 1999). 

So, we have fitted our model to the XRB introducing 
a 30% higher normalization to the data, but retaining the 
same spectral shape. According to the procedure discussed 
in Sect. 3.2, we first tried with a constant Rq obtaining 
Ro = 4.9 ±0.1 (x2 ^ i^). Again a better fit is obtained 
with a redshift dependence of R. The fit to R{z) yields 
ai = 2.8±0.2, aa = -1.5±0.1 (x^ = R{z) is shown 

in Fig. 5, while the corresponding fit is shown in Fig. 6. 
With respect to the fit with the old normalization, there 
is a more pronounced increase of the fraction of AGN2 for 
z between 0.5 and 2. 

The improvement is significant at the 98% confidence 
level, according to the F-test. It is worth noting than now 
ASCA data are well reproduced, even if not used in the fit. 
This supports the higher normalization hypothesis. Both 

are significantly higher than the fits obtained with the 
old normalization, mainly due to a clear deficiency of the 
model below ~5 keV (where, however, there may be a 
contribution from other classes of sources). In this range 
the contribution of the AGNl, which cannot change in 
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Fig. 6. Same as in Fig. 3, but introducing a 30% higher nor- 
mahzation to the Gruber (1992) data. 



our model, is most relevant, and it was already deficient 
with the lower normalization. At higher energies, instead, 
the data are well fitted thanks to the increased number of 
AGN2. However, we remind that a disagreement between 
the model and the data starts to emerge for E > 7Q keV. 
This is unavoidable, as mentioned in Sect. 3.1, when the 
cut-off energy value is fixed in order to have the XRB 
energy density peak at ^30-40 kcV. 

The most interesting result is shown in Fig. 7, illus- 
trating soft and hard counts. The higher normalization 
enables to reproduce both. Even if the normalization prob- 
lem must still be considered an open issue, it cannot help 
being noted that with the high normalization a global so- 
lution is found. 



rj^ 100 
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A 
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4. Discussion and conclusions 

In the previous section we have discussed the possibility 
of a dependence on redshift of the AGN2/AGN1 ratio. In 
particular, this ratio remains constant or slightly increases 
up to z ^2, and then decreases. This would suggest that 
AGN2 are a later evolutionary stage of the AGN phe- 
nomenon, a possibility worth to be explored theoretically. 

Another possibility is that the decreasing fraction of 
AGN2 for z ^ 2 is only apparent, and that in reality there 
is an increase of the fraction of sources with Nh ~ 10^^ 
cm~^, i.e. completely hidden at all X-ray energies. This 
could be linked with the star formation rate history, which 
is observed to increase with the redshift up to Z'~2, and 
than stays constant (Madau et al. 1996; Rowan- Robinson 



Fig. 7. Same as in Fig.5, but introducing a 30% higher nor- 
malization to the Gruber (1992) data. 

1999). A high star formation rate would imply a large 
amount of dust and gas, and then a large absorption. 

A different approach in fitting the XRB consists in a 
luminosity dependent number ratio R{L), as has already 
been done by Gilh et al. (1999b). 

The direct way to discriminate between different evo- 
lutionary models is to study the AGN2 XLF, a task within 
the capabilities of the new generation X-ray missions 
{Chandra and XMM). 

In Tab.l and Tab.2 we report the AGNl and AGN2 
densities for different flux limits corresponding to the R{z) 
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r iLix 






A.GN2 perc6ntcig6 


erg CTTi~^s~^ 


dea~''' 




% 




162.8 


242.1 


60 


5 X lO"^'-' 


42.6 


85.7 


67 




18.3 


48.3 


73 


5 X 10-" 


1.6 


6.4 


80 



Table 1. AGN2 percentage prediction as a function of the 
sampling flux in the 5-10 keV band. The second and third 
column indicate the number density of AGNl and AGN2. 



Flux 


AGNl 


AGN2 


AGN2 percentage 


-■2 -i 

erg cm s 


deg'" 


deg^'-' 


% 


10"'" 


162.8 


299.3 


65 


5 X lO"^'^ 


42.6 


115.2 


73 


10-" 


18.3 


69.2 


79 


5 X 10"" 


1.6 


11.3 


88 



Table 2. Same as in tab.l, but introducing a 30% higher nor- 
malization to the Gruber (1992) data. 



models without and with the inclusion of the 30% increase 
in the normalization of the XRB. The effect of R{z) shows 
up in an AGN2 percentage decreasing at lower fluxes, a 
consequence of the sampling at higher redshifts where the 
AGN2/AGN1 number ratio decreases. It is worth noting 
that 10^^"'' erg cm^^ s^^ and 5 x lO"'^^ erg cm^^ s^^ arc 
the flux limits expected for the Deep observations of the 
Lockman Hole and the Hubble Deep Field scheduled for 
Chandra and XMM, respectively. 
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